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Quantum point contact on graphite surface 
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Department of Physics, Bilkent University, Bilkent 06533, Ankara, Turkey. 

The conductance through a quantum point contact created by a sharp and hard metal tip on 
the graphite surface has features which to our knowledge have not been encountered so far in 
metal contacts or in nanowires. In this paper we first investigate these features which emerge 
from the strongly directional bonding and electronic structure of graphite, and provide a theoretical 
understanding for the electronic conduction through quantum point contacts. Our study involves the 
molecular-dynamics simulations to reveal the variation of interlayer distances and atomic structure 
at the proximity of the contact that evolves by the tip pressing toward the surface. The effects 
of the elastic deformation on the electronic structure, state density at the Fermi level, and crystal 
potential are analyzed by performing self-consistent-field pseudopotential calculations within the 
local-density approximation. It is found that the metallicity of graphite increases under the uniaxial 
compressive strain perpendicular to the basal plane. The quantum point contact is modeled by a 
constriction with a realistic potential. The conductance is calculated by representing the current 
transporting states in Laue representation, and the variation of conductance with the evolution of 
contact is explained by taking the characteristic features of graphite into account. It is shown that 
the sequential puncturing of the layers characterizes the conductance. 
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I. INTRODUCTION 

Graphite was a prototype sample in early experiments 1 
aiming at atomic resolution in scanning tunneling mi- 
croscopy (STM), since atomically flat and perfect sur- 
faces over thousands of angstroms can be achieved in 
air. Despite this convenient situation some STM data 
obtained from the graphite surface were rather puz- 
zling, and have been the focus of interest. For exam- 
ple, that only three alternating atoms, out of six atoms 
of hexagons in a honeycomb structure can be imaged 
demonstrated the crucial role of the sample electronic 
structure in STM. 2-4 The giant corrugation 5 recorded at 
a small tip-sample distance came out rather unexpect- 
edly, and has pointed out the importance of tip-sample 
interaction effects 6-8 in STM. Atomically resolved images 
of graphite with corrugation up to 24 A were obtained 
in air while the tip progressively displaced over 100 A . 
It was argued that the contamination under the tip dis- 
tributes the perpendicular tip force over a large area, 
and that the elastic deformation of the tip amplifies the 
corrugation. 5:9 Moreover, the conductance of the atomic 
size contact and its variation with the perpendicular dis- 
placement s of the tip have been rather confusing, 10 and 
yet cannot be reconciled with what has been learned from 
metallic contacts. 11-14 

Experiments 11 on a metal sample (such as Ag) showed 
that the conductance versus the perpendicular displace- 
ment of the tip, G[s), is initially very small, but increases 
exponentially because of tunneling through the vacuum 
gap between the tip and sample surface. Usually, it first 
saturates at G < 2e 2 n/h upon the onset of contact, 6 then 
increases discontinuously with the growth of the contact. 
Earlier, based on the first calculation of conductance of 
a three-dimensional, constriction the abrupt changes of 



G(s) were attributed to the discontinuous change of the 
cross section of contact ^4 c (s). 14 Recent experiments 15 ' 16 
achieving simultaneous measurements of the perpendic- 
ular tip force F z (s) and conductivity G(s), as well as ex- 
tensive molecular dynamics simulations, 17-20 have con- 
firmed this argument. The question of whether the con- 
ductance in an atomic size contact or in a nanowire (pro- 
duced by retracting the tip from an extensive indenta- 
tion) is quantized is a subject of current interest. 21 

The contact is usually set by a single atom at the apex, 
but it grows by additional tip atoms engaging in contact 
with the metal sample. 18 ' 19 Initially, the contact radius 
R c is small, and ~ 2 — 5 A; it is in the range of the Fermi 
wavelength Xp of metals. For alkali metals, even metals 
having a half- filled s band (Au, Cu, etc.), the concepts 
and scales adopted from a free-electron gas (FEG), such 
as spherical Fermi surface, Xf, Ep, etc., can be used 
to characterize the contact and its conductance. Under 
these circumstances, once R c ~ \p, the level spacing of 
electrons (transversally confined to the contact) can be in 
the eV range; such a quantization of electronic motion in 
the constriction is reflected by ballistic electron transport 
through the contact, even at room temperature. In devel- 
oping the theory of conductance through a metal contact, 
FEG systems with metallic densities were usually postu- 
lated. It was also implicitly assumed that the tip and 
sample have the same metal. Under these assumptions 
the conductance can be calculated by representing the 
contact by a finite (but short) constriction in which the 
potential is uniform. For such a model, even Sharvin's 
formula 22 Gs — (2e 2 / 'h)(TrR c / 'Xf) 2 is able to describe 
the overall behavior of conduction, except for features 
arising from detailed atomic structure. That such a de- 
scription of a metal contact cannot always be valid, and 
hence the nature of atoms at the contact may be crucial, 



1 



(a) 




, Ri= ax 
R 2 =a(-ix + f y) 
R 3 = c z 



a (^ x+ 27fy) 



T P'=75 a 5' + 7 cz 



a = 2.456 A (room 

c = 6.696 A temperature) 



= 2.4589 A 



(zero 



c = 6.674 A temperature) 




R . f 


,..---"A 






u 





«.=¥(*+ir) 

,, 271 

G 3 = — z 



of conductance, involve several features different from a 
metal contact. So far the tunneling between a metal tip 
and a graphite surface in STM has been treated only 
qualitatively, yet no quantitative study has been pro- 
vided. 

This paper presents a thorough analysis 25 of the for- 
mation and growth of the contact on a (0001) surface 
(or basal plane) of graphite, and provides a theoretical 
understanding of electron transport through this contact 
by taking into account characteristic features which do 
not exist in usual metal contacts. To this end, the defor- 
mation of layers and the evolution of an atomic struc- 
ture with the tip pressing against the surface is cal- 
culated by using classical molecular dynamics with an 
empirical potential 26 derived from Tersoff's many-body 
potential. 27 In this work, this potential is further elab- 
orated upon to work better under high uniaxial strain. 
The effect of the tip induced deformation is revealed by 
the self-consistent-field (SCF) pseudopotcntial calcula- 
tion of the total energy, electronic band structure, and 
total density of states as functions of the lattice param- 
eter c. The results of atomic simulations and ab initio 
calculation of the electronic structure and crystal poten- 
tial are combined in a model constriction, whereby the 
variation of conductance is analyzed by calculating the 
current within the linear response theory. 



II. ATOMISTIC SIMULATIONS 



FIG. 1. (a) The atomic arrangement and the unit cell of 
the Bernal graphite. The equilibrium lattice parameters at 
zero temperature and at room temperature are taken from 
Refs. 28 and 29, respectively, (b) The first Brillouin zone and 
its irreducible wedge with symmetry points and directions. 

becomes evident in the ferromagnetic nanowires. 23,24 In 
fact, the self-consistency of the potential, the mismatch of 
Fermi surfaces between electrodes (tip and sample) that 
are made of markedly different metals, and conservation 
of momentum have not been treated thoroughly yet. 

The situation in graphite is even more complex and 
rather different from metal contacts. From the electronic 
structure point of view, graphite is a semimetal with 
complex and narrow Fermi surface along the (HKH) 
edges of the hexagonal Brillouin zone (BZ). Energy bands 
contributing to the Fermi surface are rather flat, and 
the free-carrier density is rather low. From the atomic 
structure point of view, graphite exhibits a strongly di- 
rectional bonding structure that leads to atomic planes 
(basal planes) with wide interlayer spacings. Owing to 
the weak interlayer coupling the crystal potential at the 
interlayer region is only ~ 2eV below the Fermi level. As 
a result the transport and mechanical properties (elas- 
tic stiffness tensor, phonon spectrum) exhibit a strongly 
directional behavior. The evolution of atomic structure 
during the growth of contact, and the related variation 



Each carbon atom with its three sp 2 hybrid orbitals 
is attached to the three nearest-neighbor C atoms in 
the same plane. This way C atoms are arranged in a 
honeycomb structure in the xy plane and form the in- 
dividual (0001) basal plane, or graphene. The graphite 
crystal forms by stacking of graphenes along the z direc- 
tion (or [0001] direction) with an equilibrium separation 
d = 3.348 A at room temperature. In the normal stacking 
sequence (of Bernal graphite) one plane is shifted relative 
to the adjacent plane, so that three alternating (a) atoms 
of a hexagon directly face three atoms (a') of the adjacent 
graphenes. Accordingly, remaining three (/?) atoms face 
the centers of the hexagons (i.e. the H sites) in the ad- 
jacent graphenes. This causes ABABAB stacking. The 
atomic arrangements, lattice parameters, and the first BZ 
with symmetry points are described in Fig. 1. In rhom- 
bohedral graphite, the layer stacking is ABCABC. The 
AAA sequence of graphenes occurs in the first stage in- 
tercalates. Owing to the strong bonding combination, of 
sp 2 orbitals between two nearest-neighbor C atoms, the 
C-C distance of a hexagon is ~ 1.418 A; it is even shorter 
than the C-C distance in diamond structure. The cohe- 
sive energy is ~ 7 eV/atom 30 and the strength is rather 
high within the graphene. In this respect the graphene 
is an essential unit of graphite. On the other hand, the 
interplanar interaction and resulting binding energy is 
weak and occurs through the small overlap < p z \H \p z > > 
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between the atoms of the adjacent layers and partly 
through the long range van der Waals interaction. Pseu- 
dopotential calculations 4 ' 31 predicts total-energy differ- 
ences smaller than 5 meV/atom among three types of 
layer sequences (ABAB..., ABCABC..., AAAA...) The 
interlayer binding interaction (the experimental exfolia- 
tion energy) is only 22.8 meV/atom. 32 The present and 
previous 33 pseudopotential calculations yield very close 
energies for interlayer binding. In Sec. Ill, we show how 
these interactions are reflected to the electronic proper- 
ties. 

The empirical potential used in our molecular dynam- 
ics simulations starts from the form extended to multi- 
layer graphite inRef. 26 by combining three distinct po- 
tentials: 

ViRij) = [V T {R t3 ) + V G {R t3 )]F{R l3 ) 

+ V R (R ij )[l-F(R ij )}. (1) 

Here Rij is the distance between the ith and jth atoms. 
The total energy of atomic system for a given configura- 
tion is then expressed by Et — \ J2i^j The main 
contribution to V(Rij) is Tersoff's potential 27 Vr(Rij), 
that yields a good description of the bonding in diamond 
and graphene. Owing to the relatively shorter range, the 
weak interaction between layers are not included in Vt- 
The potential Vg is introduced to include interlayer inter- 
action, and a strong repulsive potential Vr prevents the 
layers from collapsing. The Fermi function F(Rij) pro- 
vides a smooth transition from the many-body combina- 
tion Vr(Rij)+VG(Rij) to the repulsive potential Vn(Rij). 
The above description is similar to that of the pair po- 
tential, however, there is an important difference due to 
the implicit many-body interactions in Vt and Vg- Ter- 
soff went beyond the conventional two- and three-body 
potentials in transferability and accuracy by introduc- 
ing a new scheme in view of the quantum-mechanical 
arguments brought about by the universal binding en- 
ergy (Rydberg) curve of Rose and co-workers. 34 Further- 
more, following Abell, 35 Tersoff incorporated the bond 
order as depending on the local atomic environment in 
the empirical potential. To avoid instabilities under ex- 
cessive uniaxial strain, we made the following improve- 
ment for the potential given in Eq. (1). Owing to the 
diamond-to-graphite transition data in its construction, 
the potential normally "chooses" the tetrahedral local 
environment rather than the graphitic coordination un- 
der the influence of high pressure. However, in the case 
of uniaxial compression, there appears to be no inter- 
layer interaction if the interlayer spacing is between 2.46 
and 2.87 A, because Vr(R%j) and Vc(Rij) are identically 
zero when R i3 > 2.46 A and R i3 < 2.87 A, respec- 
tively; and the range of Vr is even shorter than that of 
Vt- For this reason, the potential results in a flat re- 
gion between c = 4.92 and 5.74 A, as seen in the inset of 
Fig. 2(a), and the slope gained by the Tersoff potential 
remains too small in comparison with the experimental 
curve. 36 We made a linear interpolation of Vt and Vq 
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FIG. 2. (a) Comparison of the experimental total-energy 
variation with the improved potential. In the inset a simi- 
lar comparison is made with the potential given in Ref. 26. 
(b) The comparison of the transferability of the interatomic 
carbon potential with ab initio data. The values on the hor- 
izontal axis represents dimer (1), graphite (3), diamond (4), 
simple cubic (6), body-centered cubic (8), and face-centered 
cubic (12) structures. 

from Rij = 1.8 A to R i3 = 2.87 A by keeping the 
parameters 26 ' 27 unchanged. This way the range of V G 
is extended toward that of Vt, and hence a weak bar- 
rier is added to V(Rij) within the interpolation interval. 
By this improvement the interlayer interaction curve be- 
comes closer to that obtained by fitting to the experiment 
as shown in Fig. 2(a). The transferability of the poten- 
tial is not destroyed, and is even slightly improved for 
high-coordinated structures when compared to ab initio 
data, 37 as seen in Fig. 2(b). 

The formation of an atomic size contact on the 
graphite (0001) surface is simulated by using the classical 
molecular-dynamics method with the empirical potential 
explained above. In order to study the effect of the size 
on the results of simulations the tip-sample system is 
represented by three models of different size. The first 
one comprises six layers, and the total number of sample 
atoms is 2016. The second (third) model is comprises 
eight layers, corresponding to 3584 (4608) total number 
of sample atoms. The positions of the atoms in the last 
two layers are kept fixed. The hard, sharp metal tip 
is represented by a robust diamond tip; it comprises 13 
(111) diamond planes, and contains 167 carbon atoms. 
The apex of the tip has a single atom, and the follow- 
ing layers contain 3, 3, 6, 6, 10, 10, 15, 15, 21, 21, 28, 
and 28 atoms. The z axis is taken perpendicular to the 
(0001) surface and the periodic boundary conditions are 
imposed in the xy plane. The temperature is rescaled to 
2 K at every two steps to avoid possible divergences in 
the kinetic energy of moving atoms. The time step t is 
taken to be 10~ 16 s, in accordance with high-order Gear 
algorithm. 38 Initially the sample is equilibrated in <~ 500 
relaxation steps before the tip starts to be pushed down 
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FIG. 3. (a) Variation of the layer heights with the tip dis- 
placement s. (b) corresponding force variation. The tip is 
positioned above (3 site, (c) and (d) are the same for the H 
site. The dashed curves indicate the layer heights averaged 
over the heights of the atoms in the same layer of graphite. 
The solid curves correspond to the averaging at close proxim- 
ity of the tip. The thick line is the position of the apex of the 
tip. (e) Snapshots of the evolution of contact at the (3 site, 
(f) The same for the H site. 

from the height h = 2.5 A. As h decreases, the number 
of initial relaxation steps increases; the equilibration is 
terminated when the fluctuations in the total energy is 
settled down. After the initial relaxation stage, the tip 
is pushed at a rate of 1 x 1CT 4 A per time step for 500 
steps, and then the system is relaxed during the following 
500 steps. Even if the velocity of the tip is faster than 
the experiment, it is small enough to allow the system 
to be equilibrated between successive instabilities if any 
occur. That the average speed of the tip, v — 50 m/s is 
appropriate, and the results of simulation are converged 
are tested by performing calculations with different push- 
ing speeds and by analyzing the variation of temperature 
and potential energy of the whole system. 

To examine the formation and evolution of the contact 
the tip is located at various special positions (a, (3, and H 
site) on the graphite surface, and it is pushed toward the 
surface (along the z axis) by keeping its (x, y) position 
fixed. The effect of the tip displacement and resulting 
perpendicular tip-sample force are displayed in Fig. 3 for 
different lateral positions of the tip. In the top panels 
are illustrated the averaged change of the layer heights 
with the tip displacement. For the surface and subsur- 
face layers the average change of heights only at a close 
proximity to the tip are also shown. The bottom panels 
show the force variation F z (s). At the beginning, h = 2.5 
A; the tip-sample force is attractive in both cases. How- 
ever, when s = 1 A, for the (3 site the force enters into 
the repulsive range and increases with increasing s. The 
repulsive range, in which the layers are compressed and 
hence the interlayer distance are decreased locally, last 
until s ~ 4 A [see Fig. 3(b)]. At that point F z (s) drops 
suddenly but becomes attractive due to the puncture of 
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FIG. 4. Variation of the tip-surface separation (calculated 
at the close proximity of the tip) with the tip displacement 

the first graphite layer under high local pressure. 39 The 
attractive force, that originates mainly from the attrac- 
tion between the atoms at the side of the tip and those 
on the surface, coexists in the course of the push, but 
becomes dominant only after the release of the repulsive 
force at the apex following the puncture. The strength 
of the attractive interaction of the metal tip cannot be 
as strong as that of the diamond tip. Note that the de- 
crease of the interlayer distance could be much smaller 
than s and so the puncture would have occurred for rel- 
atively larger s if the graphite slab under study had in- 
cluded more layers. For 4 A < s < 7 A, F z (s) stays in 
the attractive range but changes to the repulsive range 
for s > 8 A, where the strong repulsive forces between 
the apex and second layer start to dominate the exist- 
ing attractive forces. This behavior continues periodi- 
cally. We note that the force variation for the H site in 
Fig. 3(d) takes place in reverse order relative to that of 
the f3 site. This is due to the fact that in the Bernal 
graphite the f3 and H sites of graphite planes occur alter- 
nately along the z axis. Therefore, when the tip starts at 
the (3 site of the first layer, it faces the H site of the second 
layer following the puncture of the first one. The impor- 
tant results obtained from these atomistic simulations are 
summarized: With increasing s, the tip first presses the 
layers. This causes the interlayer distance to decrease 
locally. Once F z (s) reaches a threshold value, the first 
layer is punctured, whereby the repulsive force is released 
and the graphite layers under compression are relaxed to 
maintain their interlayer distance, temporarily. However, 
since s continues to increase the compression of graphite 
starts again and the same events repeat quasiperiodically, 
i.e. compression of layers, puncture of the layer under 
the tip, and relaxation. The evolution of the contact is 
described by the snapshots of the atomic structure ob- 
tained from simulations in Figs. 3(e) and 3(f). Note that 
the above behavior is different from the contact formed 
on the metal surfaces, where the quasiperiodic cycles, i.e. 
compression-relaxation followed by the puncture of layer 
is absent; the contact grows by the implementation of 
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atoms to the contact area. 17-19 While the puncture of the 
atomic plane occurs through the breaking of the bonds 
at and H sites, the puncture initiate the formation of 
flakes if the contact is set at the a site. 

Another interesting feature that is absent in metal con- 
tacts is the jump to contact of the graphite surface to- 
wards the tip. As the tip approaches the sample the 
surface atoms at close proximity of the tip first move 
toward the tip and then maintain the separation h s ap- > 
proximately unaltered for a significant displacement of ^ 
the tip even if the interlayer distance decreases, (see g 
Fig. 4) . This situation lasts until the plastic deformation £ 
sets in. The attraction of the atoms under the tip result- 
ing in the jump to contact depends on the position of the £ 
tip (a, 0, and H sites) and its orientation relative to the 
honeycomb structure. 



III. ELECTRONIC STRUCTURE 
CALCULATIONS 

Effects of the tip-induced deformation on the elec- 
tronic energy structure and electronic potential of the 
contact are investigated by using SCF pseudopoten- 
tial calculations in momentum space within the local- 
density approximation and by using the Ceperley- Alder 
exchange-correlation potential. 40 The ionic potential of 
carbon is represented by nonlocal, norm-conserving 
pseudopotential. 41 The tip-induced deformation is repre- 
sented by a uniaxial strain and total energy, band struc- 
ture E(k), local density of states p(E,r), total density 
of states D(E), and crystal potential V c (r) are calcu- 
lated for different interlayer distances. The kinetic en- 
ergy cutoff (h 2 /2m)\k + G| 2 is taken to be 37 Ry that 
corresponds to approximately 900 plane waves for equi- 
librium structure. The irreducible wedge of the BZ is 
sampled by uniformly distributed 48 k points. The con- 
vergences relative to the plane wave basis set and the 
k-point sampling in BZ are tested by repeating the cal- 
culations with 33 Ry < |^|k + G| 2 < 45 Ry and with 
uniform mesh points 216, 360, 432. The change in the 
total energy is found to be smaller than 0.5% in each 
case. Taking into account the accuracy achieved by stan- 
dard LDA calculations 42 , the kinetic-energy cutoff and 
k-point sampling used in the present calculations are 
found to be appropriate to reveal the effect of deforma- 
tion on the electronic transport. The SCF cycles are iter- 
ated until rms deviation of the potential is smaller than 
10~ 7 Ry. The calculated band structure of bulk graphite 
(with equilibrium structure) is shown in Fig. 5. It is in 
general agreement with previous calculations that used 
norm-conserving, 43 ' 44 soft-transferable, 33 and ultrasoft 31 
pseudopotentials, and also full-potential all-electron cal- 
culations with linearized augmented plane wave, 45 lin- 
ear muffin-tin orbitals, 46 and linear combination of Gaus- 
sian orbitals. 47 The band energies are also in good agree- 
ment with experimental data (angle-resolved photoemis- 



18 
14 
10 
6 



f-: 



J 3 *- si. 
• • • • 



2 \j£$- 



-22 



--3^ ) •<&•: — 



"•ST 



-10 



-18 



°2i 

a 



K M 



K H 



L HH KA r L M 



FIG. 5. The energy-band structure of graphite calculated 
by the SCF-pseudopotential method. The equilibrium lattice 
parameters and corresponding Brillouin zone are illustrated 
in Fig. 1. The zero of the energy is taken at the Fermi level. 



sion spectroscopy, 48 infrared reflectance spectra, 49 ' 50 the 
angle-integrated photoemission spectroscopy, 51 and an- 
gle resolved inverse photoemission spectroscopy 52 ), ex- 
cept the a-iv and o% v bands at T-point. As expected 
the LDA results underestimate the band gap. Differ- 
ences among LDA calculations arise due to different k- 
point sampling and energy cutoff used in the calculations. 
However, the splittings at the K point, which are relevant 
for the present study are among the best LDA results that 
reproduce the experimental data. 

The energy band structure of graphene is essential for 
bands of bulk graphite. Each graphene has four filled 
bands; three bands are due to the a bonds of the lateral 
sp 2 hybrid orbitals and one band is due to the dangling 
p z orbital forming delocalizcd bonds. At the T point 
these bands are labeled by <j\ v , -K\ v , a 2v , and <7 3v in or- 
der of increasing energy The empty bands are labeled 
by al, £72, 03, and 7rJ\ The band structure of graphite, 
that has four carbon atoms in the primitive unit cell (or 
two periodically repeating graphenes) can be viewed as 
if it consists of two superimposed graphene bands which 
are slightly shifted (split) owing to the weak interlayer 
coupling. The dispersion of these bands are large (5 — 10 
eV) when k is parallel to the (0001) plane owing to strong 
overlap of the orbitals in the relatively shorter lateral C- 
C distance. However, for k perpendicular to the (0001) 
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FIG. 6. The band structure of graphite is calculated for 
various values of the lattice parameter c (or interlayer distance 
d = c/2) by using the SCF-pseudopotential method. Only the 
7r bands along the relevant directions of the BZ are shown, and 
the zero of the energy is taken at the Fermi level. 



plane (or k is parallel to the z axis) the dispersion of the 
7r bands is small (2 — 3 cV); the a nv bands are almost flat. 
While the energy gap between the occupied and unoccu- 
pied bands is ~ 6 eV at the T point, it diminishes at the 
K point and along the KH direction where the bonding 
and antibonding bands join. Accordingly, graphite is a 
semimetal. Along the KH direction 7r bands originating 
from the bonding and antibonding combination of the p z 
orbitals located at a and a' atoms (i.e 7r Q and 7r*) have 
significant dispersion due to their coupling along the z 
direction, whereas the dispersion of the 7173 and 7fg bands 
along the KH direction is rather small. 

The behavior of the 7r bands along the KH direction 
where they cross the Fermi level determine the Fermi sur- 
face and hence is essential for various physical properties 
of the bulk graphite. The Fermi surface is generated by 
a k • p extension of these bands around the HKH axis 
(which is also known as the SWMc parametrization 53 ), 
or by their Fourier expansion 54 of the LDA bands. It in- 
cludes six majority electron pockets located around the 
H points, and 18 small pockets of minority electrons. The 
majority carrier concentrations, n e ~ 2.6 x 10 18 cm~ 3 
and ri/j ~ 2.2 x 10 18 cm" 3 are reported from various 
Shubnikov-de Haas and de Haas- van Alphen experiments 
measuring the period of the magneto-resistance and mag- 
netic susceptibility 30 . The experimental values are repro- 
duced by the LDA bands. 

The effect of the uniaxial strain on the electronic struc- 
ture can be explained by a simple tight-binding picture. 
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FIG. 7. The total density of states D(E) is calculated for 
various values of the lattice parameter c. The interlayer dis- 
tance is d = c/2. The density of states near Ef is highlighted 
in the inset. 



Considering only p z orbitals in the primitive unit cell lo- 
cated at a, (3, a', and (3' atoms (see Fig. 1), a and a' 
atoms have significant overlap, whereas the overlap be- 
tween (3 and (3' atoms is rather small. Consequently, the 
TT a and 7r* bands have to be dispersive, while the 7173 and 
bands stay almost flat along the KH direction. Fur- 
thermore, the dispersion of the ir a and 7r* bands along 
the KH direction increases with uniaxial compression 
since the overlap between ir a and 7r* atoms increases. 
The modification of the electronic structure under the 
uniaxial strain is illustrated by the ab initio bands shown 
in Fig. 6. In comply with the above arguments, the ir a 
and 7T* bands are strongly affected from the compres- 
sive strain. The dispersion of the doubly degenerate 7173 
bands in Fig. 6 is small but essential in determining the 
size of the electron pockets of the Fermi surface. Note 
that already with compression the height of the BZ (i.e. 
HKH distance in the k space) increases. Furthermore, 
the overlap of the it bands along TK direction increases 
with increased compressive strain. For c ~ 4.67 A the 
7173 band rises and touches the Fermi level along the TM 
direction. 

It becomes clear from the above discussion that the 
metallicity of the bulk graphite increases with increasing 
uniaxial compressive strain along the z direction. This 
situation can be revealed by examining the total density 
of states at the Fermi level D(Ep) for various values of 
c. Here we calculated the histogram of the state distri- 
bution by using 198 k points in the irreducible wedge of 
the BZ and we then broaden it by a 0.5 eV Gaussian 
convolution. Our results are illustrated in Fig. 7. It is 
seen that D(Ep) may increase as much as 300% upon 
~ 30% contraction of c. Such an increase at D(Ep) have 
important implications in the transport properties. 55 
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FIG. 8. Contour plots of the crystal potential on a perpen- 
dicular plane and the plane- averaged crystal potential V c {z) 
as a function of c. The contour interval is 1.5 eV in each 
graph, and the indicated contour values are the maximum 
values with respect to the Fermi level. 

Finally, we examine the variation of the crystal poten- 
tial V c (r) with a uniaxial compressive strain, i.e. with 
the interlayer spacing that is calculated selfconsistently. 
The potential counterplots are shown in a relevant per- 
pendicular plane in Fig. 8. The planar averaged crystal 
potentials (V c (z) — ^ J n V c (x,y, z)dx dy, O being the 
cross section of the unit cell at the basal plane) are also 
shown in the same figure. The crucial conclusion one can 
draw from these contour plots is that the potential be- 
tween layers are lowered (becomes more attractive) with 
decreasing c. Note that the crystal potential between lay- 
ers is high and behaves as a barrier, but it is still below 
Ep. In concluding this section, we emphasize that under 
the uniaxial compressive strain, that makes the inter- 
layer distance smaller, the Fermi surface grows, D(Ep) 
increases and V c (r) between layers is lowered and be- 
comes more attractive. Thus the metallicity of graphite 
increases. Because of the low carrier concentration and 
weak screening, such effects are expected to survive at 
the contact region. 



IV. CONDUCTANCE THROUGH AN 
ATOMIC-SIZE CONTACT 

In the electron transport between the STM tip and 
sample one normally distinguishes two regimes, i.e. tun- 
neling and ballistic. The topological mode of STM (in 
which the separation is large enough to allow a vacuum 
barrier <I>, and the bias voltage Vb is small) operates in 
the tunneling regime. Since the tunneling conductance is 
proportional to the local density of states of the sam- 
ple at the center of the tip and at the Fermi energy, 
p(Ep, Rt), 56 the tunneling current has been rather small 
in STM measurements on the graphite surface. A sig- 
nificant value for the conductance can be achieved only 
at a small tip-sample separation, where the electronic 
contact is already established. Earlier, it was found that, 
due to the tip-sample interaction, <I> collapses 6 for h ~ 3.5 
A and valence bands (normally below Ep) begin to over- 
lap with the Fermi level in regions different from the KH 
direction of the BZ. 8 Based on the ab initio electronic 
structure calculations it was also shown that the tip- 
induced states (TILS)(Ref. 6) form at the close proximity 
of the tip. Actually, TILS are nothing but the states that 
arc confined in the constriction and are precursors to the 
current transporting states. In the present work, it is also 
shown that under the uniaxial stress (which is induced 
by the tip) the valence bands, which arc normally be- 
low Ep, begin to cross the Fermi level. The Fermi-level 
crossing may induce sudden changes in physical proper- 
ties, in particular in the tunneling conduction. 57 Soler et 
al. 5 were able to obtain only G — 0.05(2e 2 //i) when the 
largest corrugation of line scans were achieved. The sig- 
nificant tip-induced deformation suggests that the barrier 
collapses at the contact. That the maximum conduc- 
tance measured is much smaller than 2e 2 /h disregards 
the opening of the first ballistic channel. Such a situation, 
in which the ballistic conductance cannot set in despite 
the collapse of the barrier, i.e., $ < Ep, was pointed out 
earlier. 58 The effective barrier $ e ff, which may be greater 
than Ep due to the size of the constriction, hinders the 
ballistic transport. 58 On the other hand, it was argued 
that the flake of a dirt between the apex of the tip and 
the graphite surface increases the resistance and amplifies 
the corrugation 9 . Nevertheless, apart from the spreading 
resistance, the conductance of a contact created by the 
sharp tip is generally smaller than 2e 2 /h. This is due to 
the small density of states D(E), available at Ep, since 

the current / oc J^ +eVB D t (E)D(E + eV B )T(E)dE. In 
this expression D t is the density of states of the tip, and 
T(E) is the transmission coefficient. The crystal poten- 
tial between graphenes is also rather shallow, only ~ 2 
eV below Ep. This may prevent a current transporting 
state (or a ballistic channel) from forming at the con- 
tact between two electrodes. As a result the transmis- 
sion coefficient T(E) becomes small. Another important 
feature that was not taken into account extensively in 
earlier works investigating the ballistic electron transport 
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FIG. 9. (a) A constriction model describing the point con- 
tact created by two Al(lll) tip on both sides of a thin graphite 
slab. The variation of the potential relative to the vacuum 
level is shown below. The potential has a cylindrical symme- 
try, (b) Variation of the conductance with the displacement 
of the tip. The dotted and dashed curves correspond to cases 
(i) and (ii) explained in the text, respectively. 4>^ is the work 
function, and r c , A c , zr, and zl are defined in the text. 

between two free-electron systems through a constriction 
is the actual topology of the Fermi surfaces. While the 
Fermi surface of a metal tip can be taken (but not always) 
as a sphere, the Fermi surface of graphite is far from be- 
ing similar to a sphere. Since the narrow Fermi surface 
occurs at the corners of the hexagonal BZ, electrons at 
Ep have a very large lateral momentum Kkx- Even if 
the lateral momentum of the incoming electron cannot 
be conserved at the entrance of the contact, it should 
gain Kk.K after it passed to the graphite, which can be 
supplied by the acoustic phonons. Therefore, the match- 
ing of the Fermi surfaces of the electrodes at both sides 
of the contact is important; the conductance measured 
at a low phonon population (or at low temperature) is 
expected to reveal interesting features. 

Here we first present a qualitative discussion of the 



conductance through a contact created by a sharp tip. 
As described in Sec. II, at a separation h s the surface 
layer of graphite is first attracted and jumps to the con- 
tact. After this stage the separation between the tip and 
surface remain approximately unchanged for a significant 
displacement of the tip, unless a plastic deformation sets 
in (see Fig. 4). Accordingly the current under constant 
Vb docs not increases significantly. As s increases, the 
tip presses the underlying graphite surface as shown in 
Fig. 3. This decreases interlayer distance, lowers the in- 
terlayer potential and increases the density of states at 
Ep (see Figs. 6-8). Even if T(E) (that also increases 
with s) is assumed to be unchanged, the current and 
hence G increase continuously with s in the elastic re- 
gion. However, as seen in Fig. 3, at certain range of s 
(where G reaches its highest value) the compressive stress 
exceeds a threshold value whereby the tip punctures the 
first layer by atomic rearrangements or by creating flakes 
at its close proximity. This relaxes the strain in the lay- 
ers. As a result, the interlayer distance and the density 
of states are retrieved. Since the tip faces now the sec- 
ond layer but with a larger separation, the conductance 
G decreases suddenly. As the tip first presses than punc- 
tures the graphite layers sequentially with increasing s, 
the conductance first increases from a low level to a high 
level and then falls again to a low level. Due to the elec- 
tron transport from the sides of the tip to the already 
punctured layers the low level value of G may increase 
with increasing number of layers. Nevertheless, G(s) is 
expected to vary between low and high conductance val- 
ues. This situation may be different if a blunt (and flat) 
metal tip presses the graphite surface and forms a contact 
with large A c . Whether the quantum ballistic transport 
occurs for large contact area depends on the matching of 
the Fermi surfaces. 

We now quantify the above discussion about G(s) by 
calculating the conductance of the contact described in 
Fig. 9. Two metal tips (Al oriented along the [111] direc- 
tion) with a single atom at the apex make point contact 
from both surfaces of the graphite (0001) plane. This 
contact can be described by a tapered constriction that 
has cylindrical symmetry and is connected to the free 
electron systems from both sides. The Al electrodes are 
taken as free-electron systems with an electron density 
and Ep appropriate for Al. The radius of the constric- 
tion through the graphite is determined from our ear- 
lier ab initio calculation dealing with Al tip and graphite 
surface 8 and present calculations. The crystal potential 
V c (r) of graphite between two adjacent layers raises to 
~ 2 eV below Ep. This creates a saddle-point effect. 6,59 
The potential of the constriction is described in Fig. 9(a). 
After the mechanical contact is established, the displace- 
ment of the tip decreases the interlayer distance c/2 lo- 
cally by applying a uniaxial compression. This, in turn, 
induces changes in the electronic structure and the crys- 
tal potential. For example, for the displacement s ~ 0.4 
A (where the puncture sets in), the saddle point poten- 
tial is lowered by 1.8 eV for a three layer slab shown in 
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Fig. 9(a). In this simple constriction model, where the 
electrodes are represented by a free-electron gas system, 
we can consider only the variation of graphite potential 
and the radius of the constriction, i.e. V c (s) and r c (s). 
Therefore, the calculation of conductance with V c (s) and 
r c (s) is expected to describe the essential features of the 
G(s) curve. 

To calculate conductivity, we used the recursion trans- 
fer matrix method proposed by Hirose and Tsukada; 60 
and it was applied in similar types of calculations in 
Ref. 61. In this approach, we replace our constriction 
shown in Fig. 9(a) by a periodically repeating supercell 
structure (60 x 60 a.u. 2 ) in the xy plane defining the con- 
tact. This allows us to express the jth scattering state in 
terms of Laue expansion: 

*j(r) = c 4k ii' r ii 5^lM*)e <G 'i' p 'i. (2) 
l 

Here k|| and G|| correspond to the wave vector in the 
BZ and the reciprocal lattice vector of the supercell in 
the xy plane, respectively. The summation over G|| is 
truncated and a cutoff energy of 30 eV was sufficient for 
convergence. Deep in the Al structure, the wave function 
*j(r) can be expressed as: 

= Eittf | ^ z | c^+S)-!, {z > Zr1 
< * 3 -(r) =e ik i-V (k|l+G ll ) - r " 

+ E 4 ^|° k : k ;' Z |c 4(k " +G j ) - r " {z<zl), 

(3) 

where zr and zl are the boundaries of the constriction. 
tij and Tij are transmission and reflection coefficients, 
corresponding to states k z and k J z , and they are deter- 
mined by using transfer-matrix method. We note that 
in this approach, the potential should be smooth and 
the constriction is divided into "small strips" along the 
propagation direction (for recursive purposes), otherwise 
results diverge. To smooth out the potential at both ends 
of Al tips, we used a Fermi function of width 3 a.u., and 
discretization width of 0.1 a.u. was applied to obtain 
convergent results. In our model, we assume that V c (s) 
varies linearly (9.3 —> 7.54 eV) as a function of s. Fur- 
thermore the mean width of the saddle-point potential 
A c also decreases linearly with s. The radius of the con- 
striction, r c (s), is also a difficult parameter to determine. 
Here we assume that r c (s) is uniform between the Al tips, 
and consider two limiting cases: (i) r c (s) = r = 2.2 A, 
i.e. the value corresponding to h s ~ 2 A at the electronic 
contact, 7 (ii) r c (s) = tq + 1.23.s, allowing a 50% expansion 
of contact area before puncture occurs. Furthermore, lin- 
ear response theory is assumed and we obtain: 

ij ij 



where i and j run over all conducting states, Tjj = 
tij x (kyki) 1 / 2 , and R y - = x (k^/k^) 1 / 2 . We have 
considered only the case where k|| = 0, since the su- 
percell area of the boundary condition is large and and 
hence the area of supercell BZ is negligible. We have 
applied both formulas in Eq. (4) for determining G, be- 
cause they serve as a good test of whether the results 
have converged. The calculated G(s) curves are shown 
in Fig. 9(b). As argued above, G(s) increases with in- 
creasing s until a plastic deformation (or puncturing of 
the graphene) sets in. Once the graphene is punctured by 
the sharp tip, the conductance decreases suddenly. For 
a sharp tip having a single atom at the apex, the radius 
of the constriction is normally small. In Fig. 9(b), G at- 
tains a value smaller than 2e 2 /h, corresponding to case 
(i) with r c = tq. For case (ii), where r c (s) is allowed 
to expand, a ballistic channel is opened before the punc- 
turing. We also note that AG/ As of actual graphite is 
expected to be much smaller than that of Fig. 8(c), since 
the same level of compressive strain in Fig. 8(c) requires 
much larger s in the multilayer graphite. For the con- 
striction described in Fig. 9(a), s ~ A c . 



V. CONCLUSIONS 

In this work, based on the atomistic simulation of the 
contact and ab initio electronic structure calculations of 
graphite under uniaxial compressive strain we showed the 
following: (i) The graphenes are punctured sequentially 
by a sharp tip pressing the surface, (ii) The tip induced 
elastic deformation between two consecutive puncture in- 
creases the density of states at the close proximity of the 
tip, and lowers the crystal potential: Hence, the metallic- 
ity of graphite increases with uniaxial compressive strain 
of basal planes, (iii) Accordingly the conductance be- 
tween two consecutive puncture increases, but falls sud- 
denly upon the onset of a new puncture. This is a be- 
havior different from the usual metal contacts. 
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